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7— I . Abstract 

> _ 

This paper is a continuation of hepth/0507224 where open topological B- models describing 
D-branes on 2-cycles of local Calabi-Yau geometries with conical singularities were studied. After 
a short review, the paper expands in particular on two aspects: the gauge fixing problem in the 
reduction to two dimensions and the quantum matrix model solutions. 
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1 Introduction 

Singular Calabi-Yau manifolds represent one of the most interesting developments in string 
compactifications. For instance, the presence of a conifold point, pQ, in a Calabi-Yau opens new 
prospects: in conjunction with fluxes and branes it may allow for warped compactifications, which 
5_i ■ in turn may create the conditions for moduli stabilization and for large hierarchies of physical 



scales. On the other hand singular Calabi-Yau compactifications with conical singularities seem 
to realize favorable conditions for low energy theory models with realistic cosmological features. 
A conifold singularity can be smoothed in two different ways, by means of a 2-sphere (resolution) 
or a 3-sphere (deformation). This leads, from a physical point of view, to a geometric transition 
that establishes a duality relation between theories defined by the two nonsingular geometries 
(gauge-gravity or open-closed string duality) , [2J E| ■ In summary, conifold singularities are at the 
crossroads of many interesting recent developments in string theory. It is therefore important 
to study theories defined on conifolds, i.e. on singular non compact Calabi-Yau threefolds, as 
well-defined and (partially) calculable models to approximate more realistic situations. 

In [31, building on previous literature, we started to elaborate on an idea that is receiving 
increasing attention: how data describing the geometry of a local Calabi-Yau can be encoded, via 
a topological field theory, in a (multi-)matrix model and how they can be efficiently calculated. 
The framework we considered was IIB string theory with spacetime filling D5-branes wrapped 
around resolving two-dimensional cycles. This geometry defines a 4D gauge theory, [51 Itjl 171 PUj. 
On the other hand one can consider the open topological B model describing strings on the 
conifold. The latter has been shown long ago by Witten to be represented by a six-dimensional 



1 



holomorphic Chern-Simons theory, When reduced to a two-dimensional cycle this theory 

can be shown to reduce to a matrix model. In particular, if one wishes to represent a wide class 
of deformations of the complex structure satisfying the Calabi-Yau condition, one must resort 
to very general multi-matrix models, [S]. In 0] we concentrated on the topological string theory 
part of the story, ^3 El El El El El? m particular on the formal aspects of the reduction 
from the six-dimensional holomorphic Chern-Simons theory to a two-dimensional field theory 
and on the analysis of the matrix model potentials originated from the Calabi-Yau complex 
structure deformations. Finally we concentrated on the subclass of matrix models represented 
by two-matrix models with bilinear coupling (for a very recent development see 116)). In this case 
the functional integral can be explicitly calculated with the method of orthogonal polynomials. 
Using old results we showed how one can find explicit solutions by solving the quantum equations 
of motion and utilizing the recursiveness granted by integrability. All the data turn out to be 
encoded in a Riemann surface, which we called quantum Riemann surface in order to distinguish 
it from the Riemann surface of the standard saddle point approach. 

In this paper we would like to return to some topics that were only partially developed in 
In particular, in section 2, after a concise review of the reduction to from CS theory to matrix 
models, we explain in detail the gauge fixing problem in this process. Subsequently we return 
to the problem of solving two-matrix models with bilinear couplings by means of the orthogonal 
polynomials method via the solution of the quantum equation of motion. After a short review of 
the method in section 3, our main purpose in section 4 is to clarify the similarities and differences 
with other methods, in particular with the semiclassical saddle point method. We discuss at 
length the result of [3] that the quantum equations of motion admit in general more vacua than 
the saddle point method. We interpret these additional solutions as 'quantum' cycles that have 
no classical analog. Finally in section 5 we give a simple, explicit example of topological open 
string expansion. 



2 Reduction to the brane 

In this section we summarize the reduction of the topological open string field theory (B model) to 
a holomorphic 2-cycle in a local Calabi-Yau threefold [I]. Let us consider a holomorphic C 2 — fibre 
bundle X —* S, where £ is a Riemann surface. The space X is obtained as a deformation of the 
complex structure of the total space of a rank-2 vector bundle V on T,. Given an atlas {U a } on 
S, the transition functions for X can be written 



Z(a) = f(a/3)(z(/3)) 



»,j = l,2 



where /(«/?) are the transition functions on the base, Mj/m the transition functions of the vector 

bundle V and ^ 3 r a p\ are the deformation terms, holomorphic on the intersections (U a fl Up) x C 2 . 

The Calabi-Yau condition on the space X, i.e. the existence of a nowhere vanishing holomor- 
phic top-form £1 = dz A dw 1 A dw 2 , puts conditions on the vector bundle and on the deformation 
terms. The determinant of the vector bundle has to be equal to the canonical line bundle on £ 
and for the transition functions this means detM^) x ft a p\ = 1- For the deformation terms 

we have det (1 + d*&) = 1, where (1 + d^>) l j = 6) + dj^f 1 . The solution of this condition can be 
given in terms of a set of potential functions X( a p^, the geometric potential, which generates the 
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deformation via 



where we define the singular coordinates wi r a p) = + ^\a/3) i^iP)) • 

The topological open B-model on X can be obtained from open string field theory and reduces 
jllj to the holomorphic Chern-Simons (hCS) theory on X for a (0,l)-form connection on a U(N) 
bundle E, where N is the number of space-filling B-branes. We will restrict ourselves to the case 
in which E is trivial. The action of hCS is 

S(A) = — [ C, C = tlATr ( -AAdA+-AAAAA) (2.1) 
9s Jx V 2 3 / 

where A £ T^ 0,1 \X). The dynamics of B-branes on a 2-cycle ScX can be described by reducing 
the open string field theory from the space X to the B-brane world- volume £. 

To obtain the reduced action for the linear geometry (vl/j = 0), first we split the form A into 
horizontal and vertical components using a reference connection T on the vector bundle, then we 
impose the independence of the fields on the vertical directions and finally we "integrate along 
the fiber" using a bilinear structure K on the bundle. If the connection V is the generalized 
Chern connection for the bilinear structure, then the result is independent of the particular 
(r, K) chosen. 

Let us define Az = A s — A^F^w 3 and Ai = Aj, where T is a reference connection and impose 

that the components (Ag,Aj) are independent on the coordinates along C 2 , obtaining for the 
Lagrangian 

L = -U A Tr [A-DzAj + A-T^Aj^ dw J Adz A dw~ j (2.2) 

where Dz is the covariant derivative w.r.t. the gauge structure. 

Now let us consider a bilinear structure K in V, i.e. a local section K S T(V ® V), the 
components K 13 being an invertible complex matrix at any point. The "integration along the 
fiber" is realized contracting the hCS (3,3)-form Lagrangian by the two bi-vector fields k = 



Cred = ipAkL = -dzdz(det K)e ij Tr AjDgAj + AjT^A^ . (2.3) 



Defining the field components tp 1 = iy^A £ V, where V 1 = K K '§-=7, one gets 
Cred = ^dzdzTr 



eirfDrf + (det K)<p m W ( K^ t d- z K n3 + K^K nl Y% ) (2.4) 



where Kjj are the components of the inverse bilinear structure, that is K^jK jl = 6k. In order 



to have a result which is independent of the trivialization, just set the reference connection to 
be the generalized Chern connection of the bilinear structure K, that is = KjidzK lk . The 
action for the reduced theory is given by 

Sred = — [ Cred = TT" / dzdzTr [cijCf 1 Dglf>] . 

9s Jt, l 9s Jt, 

In the rational case S ~ P 1 with non vanishing deformation terms ^j, X is a deformation of 
a vector bundle V — Opi (n) © Opi (— n — 2) for some n. 
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Let us start with the reduction in the Abelian U (1) case. In this case the cubic term in the hCS 
Lagrangian is absent and the reduction is almost straightforward. In the singular coordinates 
(ip 1 ,^ 2 ) we obtain that 



C 



red 



2 6ij 



ip l dz^f 3 dzdz 



(2.5) 



in both charts of the standard atlas {Un, Us} of P . The potential term X gives the deformation 
of the action due to the deformation of the complex structure. Passing to the non singular 
coordinates (4>i,4>2), one obtains 



'red 



1 

2g s 



dzdzeijdfdzd) 3 + 



dz 
2ttI 



X(z,, 



(2.6) 



where § is a contour integral along the equator. Therefore, the reduced theory gives a b-c (J3—^) 
system on the two hemispheres with a junction interaction along the equator. 

The non- Abelian case is a bit more complicated than the Abelian one because of the tensoring 
with the (trivial) gauge bundle. This promotes the vector bundle sections to matrices and 
therefore it is not immediate how to unambiguously define the potential function X in the 
general case. The easiest way to avoid matrix ordering prescriptions is to restrict to the case 
in which X(z,ui) does not depend, say, on to 2 . Defining B := uj 1 ^ 2 + X, one obtains \E' 1 = 0, 
^ 2 = d,„iB and the reduced action 



S=S, 



1 



red 



-Tr(<f) 2 D z -(f> 1 )dzdz + jltBiz,^)^ 



(2.7) 



2.1 Gauge fixing 

In order to calculate its partition function, let us now discuss the gauge fixing of the theory. The 
following discussion is a refinement of the derivation given in [5]. Our starting action is Q2.7JI and 
we follow the standard BRST quantization (see for example 8 ). 
The BRST invariance in the minimal sector is 

sA s = -(Dc) 2 , s4> % = [c,4> 1 ], sc = -[c,c] 

while we add a further non minimal one to implement the gauge fixing with 

sc = b, sb = 0. 

The gauge fixed action is obtained by adding to S a gauge fixing term 

S af = S + s*, where * = — / Ticd z A z 

9s J pi 

which implements a holomorphic version of the Lorentz gauge. Actually we have 

s ^ = - [ Tr [bd z A- z - d z c(Dc) z ] . 
9s J pi 

Our partition function is then the functional integral 

Z B = fv 14>\A- Z ,c,c,b) e~ s of 
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The calculation can proceed as follows. Let us first integrate along the gauge connection A 2 
which enters linearly the gauge fixed action and find 

Z B = Jv [cft\c,c,b] e ~ii[-Ip^^-^^f^M 1 )]s{ dzb+ [ dz c,c] + [<P\0 2 ]} 



Now we integrate along the field b. By solving the constraint we obtain 



Z B = J V[(ft\c,c]e 



Ml 1 



det'a 



where det' is the relevant functional determinant with the exclusion of the zero modes. Then we 
integrate along the (c, c) ghosts and get 

Z B = [v U*] e~T s [- ^ T*B{z&)) det'd 2 d s 

J det'd z 

Finally, since the geometric potential B does not depend on (ft 2 , we can also integrate along this 
variable and obtain 

The delta function constrains the field (ft 1 to span the c^-zero modes and once it is solved it 
produces a further (det'<9g) multiplicative term that cancel the other determinants. Therefore, 
all in all we get 

z B = I 

JKerdz 

Lastly we can expand (ft 1 = X^I=o^^ along the basis £i(z) ~ z % of Ker<9 5 with N x N matrix 
coefficients Xj. Finally we find the multi-matrix integral 



Z B = I f[dX ie -^ Xo >~> X ^ (2i 



i=0 

where we defined 

W(X ,...,X n )= lTrB(z,J2 X i^) ( 2 - 9 ) 



This is the result of our gauge fixing procedure which covers the details needed to complete 
the derivation given in [I] and confirms the conjecture in [5]. 



3 General properties of two— matrix models. 

The second part of this paper is devoted to solving some of the matrix models introduced above, 
ea, (|2.8|) . the two-matrix models with bilinear coupling. The sequel is a short review of 4 on 
this subject containing some additional remarks and complements. The main point we insist 
on is that for these models there is the possibility to solve the quantum problem exactly. That 
is, we perform the path integral exactly and determine all the (quantum) solutions. The exact 
solvability of these two-matrix model is a well-known fact, but its consequences have not yet 
been completely explored. As we will see, not all these solutions have a classical analog and, 
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thus, they represent genuine new quantum solutions. To this purpose the method of orthogonal 
polynomials turns out to be particularly fit. This method allows one to explicitly perform the 
path integration, so that one is left with quantum equations of motion and the flow equations of 
an integrable linear systems. The latter in particular uncover the integrable nature of two-matrix 
models, which stems from the Toda lattice hierarchy underlying all of them. Our approach for 
solving two-matrix models consists in solving the quantum equations of motion and, then, using 
the recursiveness intrinsic to integrability (the flow equations), in finding explicit expressions for 
the correlators. An alternative method is based on the W constraints on the functional integral. 
We do not use it here, but one can find definitions, applications and comparisons with the other 
methods in |21l 122 j . For general references on matrix models, see the bibliography in 2£>1I3]. 

Let us start with a synthetic review of the approach based on orthogonal polynomials. The 
model of two Hermitian N x N matrices M\ and M 2 with bilinear coupling is defined by the 
partition function 

Z N (t, c) = J dM x dM 2 e trW , W = Vi + V 2 + cM x M 2 (3.1) 

with potentials 

Pa 

V a = Y j i a , r M r a a = 1,2. (3.2) 

r=l 

where p a are finite numbers. We denote by M P1 , P2 the corresponding two-matrix model. With 
reference to ea. l|2.9|) . this model descends from the geometric potential B defined by 



B(z,u) = - V 1 {uj) + V 2 [-) (3.3) 



We are interested in computing correlation functions (CF's) of the operators 

T k = trM?, a k = trM 2 , Vfc, 

For this reason we complete the above model by replacing ()3.2|) with the more general potentials 

00 

V a = J2 t <x,rK, a = 1,2 (3.4) 
r=l 

where t a>r = i Qjr for r < p a . The CF's are defined by 

Qn+m 

<r ri ... r rn <r Sl ...a Sm >= — ^-^ - — lnZ N (t, g) (3.5) 

where, in the RHS, all the t a ^ r are set equal to t a>r for r < p a and the remaining are set to zero. 
From now on we will not distinguish between t ar and t ar and use throughout only t ar . We 
hope the context will always make clear what we are referring to. 

We recall that the ordinary procedure to calculate the partition function consists of three 
steps [2Sj , [23] , j2^j : (i) one integrates out the angular part so that only the integrations over the 
eigenvalues are left; (ii) one introduces the orthogonal monic polynomials 

£n(Ai) = A™ + lower powers, rj n {\ 2 ) = \ 2 + lower powers, n = 0, 1, 2, . . . 
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which satisfy the orthogonality relations 



dX^U^e^^^+^rimi^) = K(t, c)5 nm (3.6) 

(Hi), using the orthogonality relation (|3.fij) and the properties of the Vandermonde determinants, 
one can easily calculate the partition function 

N-l 

Z N (t, c) = const N\ JJ hi (3.7) 

i=0 

whereby we see that knowing the partition function amounts to knowing the coefficients h n (t,c). 
The crucial point is that the information concerning the latter can be encoded in the flow equa- 
tions of the Toda lattice hierarchy and the quantum equations of motion. Before coming to this, 
let us introduce some convenient notations. In the sequel we will meet infinite matrices Mij with 
< i,j < oo. For any such matrix M, we define 

M = H^MH, Hij = hiS i:j , M i:j = M,i 

We represent such matrices in the lower right quadrant of the plane. They all have a band 
structure, with nonzero elements belonging to a band of lines parallel to the main descending 
diagonal. We will write M £ [m, n], if all its non-zero lines are between the m-th and the 
n-th ones, setting m = for the main diagonal. We refer to such types of band matrices as 
Jacobi matrices. Moreover M + will denote the upper triangular part of M (including the main 
diagonal), while M_ = M — M+. 

Let us come now to the quantum equations of motion. They are written as 

P°(l) + y 1 '(Q(l)) + cQ(2) = 0, cQ(l) + y 2 '(Q(2))+P°(2) = 0, (3.8) 

In these equations Q(1),Q(2),P°(1),P°(2) are infinite Jacobi matrices. They represent the 
multiplication by Ai, A2 and the derivative by the same parameters, respectively, in the basis of 
monic polynomials. Eas. (|3.8[) can be considered the quantum analog of the classical equations 
of motion. The difference with the classical equations of motion of the original matrix model is 
that, instead of the N x N matrices M\ and M2, here we have infinite Q(l) and Q(2) matrices 
together with the quantum deformation terms given by P°(l) and 7 3 °(2), respectively. From the 
coupling conditions it follows at once that 

Q(a) G [-m a ,n a ], a = 1,2 

where 

mx=P2 — l, 7712 = 1 71i = l, 71 2 =Pl — 1 

where p a> a = 1, 2 is the highest order of the potential V a (see (|3,2|0 . 
The flow equations of the Toda lattice hierarchy are 

Q(1) = [Q(1), Q k (a)-], a = 1,2 (3.9a) 



dt a ,k 
d 



Q{2) = [Q k (a) + , Q(2)} (3.9b) 

Ct a ,k 
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Finally one must use the reconstruction formula for the partition function 



a N-l 

a 



dt n 



i=0 



hxZ N (t,c) = (Q r (<x)) .., a =1,2 (3.10) 



It is evident that, by using the equations (|3.9al3.9b|) above we can express all the derivatives of 
Zn in terms of the elements of the Q matrices. For example 

d 2 ( \ 

-]nZ N (t,c) = (Q r (a)) , a = 1,2 (3.11) 



dti t idt a>r ^ / N,N-1 

and so on. We recall that the derivatives of F(N,t, c) = \nZ^{t, c) at t a>r = i a>r are nothing but 
the correlation functions of the model. 

To end this section, we collect a few formulas we will need later on. First, we will be using 
the following parametrization of the Jacobi matrices 

i 1=0 i 1=0 

where 1+ = J2i=o -^M+l anc ^ (Ei,j)k,i = ^i,k^j,l- O ne can immediately see that 

(Q + (l)) = 6 j>i+1 + aoftSij, (Q_(2)) = R(i)S jti -i (3.13) 

where R(i + 1) = As a consequence of this parametrization, ea. (|3.11|) gives in particular 

the two important relations 



d 2 d 2 
F(N,t,c) = ai (N), — - 



, F(N,t,c)= ai (N), ———F(N,t,c) = R(N) (3.14) 



To complete this summary, we should mention the W-constraints method. The latter are 
constraints on the partition function which take the form of a nice algebraic structure, see [211 122j 
for instance. They are obtained by putting together quantum equations of motion and flow 
equations. W-constraints (which are also called loop equations or Schwinger-Dyson equations) 
can be used to solve matrix models, but such a procedure is less efficient than the one used here. 

Finally we must recall that the CF's we compute are genus expanded. The genus expansion 
is strictly connected with the homogeneity properties of the CF's. The contribution pertinent 
to any genus is a homogeneous function of the couplings (and N) with respect to appropriate 
degrees assigned to all the involved quantities, see 0]. In particular we expect the partition 
function to have, in the large iV limit, the expansion 

oo 

F = ^N 2 - 29 F g (3.15) 

9=0 

where g is the genus. Such expectation, based on a path integral analysis, remains true in our 
setup due to the fact that the homogeneity properties carry over to the Toda lattice hierarchy. 
In (|3,15|) F g is interpreted as the result of summing the open string partition function at fixed g 
over all boundaries. Therefore it is a closed string quantity and the correlators obtained from it 
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are closed string correlators, i.e. correlators of the geometrically dual closed strings relative to 
the deformed geometry. From matrix models it is however possible to obtain also genuine open 
string quantities. To this end we must change expansion. We introduce the string coupling g s 
and rescale all couplings in W, so as to extract an overall factor l/g s out of them. Now the open 
string expansion is postulated to be 

oo 

F = J29s 9 ' 2+h N h F 9ih (3.16) 

h=0 

Fg^ refers to the contribution from a world-sheet of genus g with h boundaries. 

4 Solving two— matrix models 

Our procedure to solve two-matrix models consists in solving the quantum equations of motion. 
This allows us to determine the 'lattice fields' aj(n),6j(n) and R(n). Once these are known we 
can compute all the correlation functions starting from 1)3.10(1 by repeated use of eas. (|3,9al3,9b|) . 
which form the flows of the Toda lattice hierarchy. As for the free energy F(t, N, c) itself, see 
for instance ( 21 ). In the following we describe some explicit examples of this method. In 
reality we will concentrate on solving the equations of motion, since the calculus of correlators 
is of algorithmic nature and, therefore, not particularly interesting; in any case, it has already 
been illustrated in a number of examples, |22[ 1211 The equations of motion are definitely 
more interesting, because some aspects of them have not been stressed enough or ignored in the 
existing literature. Therefore our purpose is to use some (simple) examples first of all in order to 
exemplify our method and compare it with others, in particular with the saddle point method; 
second, in order to single out the novelties with respect to the existing literature. Our exposition 
will be very close to 0], but with more examples and more details. For other approaches to 
matrix models see, for instance, |26| and references therein. 

4.1 The Gaussian model 

The bi-Gaussian model A4 2 ,2 was fully solved in [22]. It is a very simple model, but we 
further simplify it by considering the decoupling limit c — > and keeping only half of couplings, 
say tk = ti ; k- This operation makes sense and leads to the Gaussian one-matrix model M. 2 - 
It allows us to make an explicit comparison of our method and the traditional one based on 
eigenvalue density and resolvent. 

The quantum equations of motion of the simplified model are, see |22| . 

ti + 2t 2 a (n) = 0, n + 2t 2 a 1 {n) = S (4.1) 

We can set t\ = 0, because the linear term in the potential can always be gotten rid of with a 
matrix redefinition. Therefore we have ao = and a\ = —n/(2t 2 ) = g s n/N. In the last equality 
we have introduced the string coupling g s = —N/(2t 2 ), thus making contact with the notation 
of ^3] where l/g s is factored out in front of the potential in the path integral. Since the latter 
reference contains a compact review of the Gaussian one-matrix model with the saddle point 
method, we refer to it for a comparison. In the large iV limit, n/N becomes the continuous 
variable x, and so, in this limit, ai(x) = g s x. Now we can easily compute all the correlators. For 
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a comparison with ^3] we need the one-point functions in genus 0. Following [2J and section 
5.5 of we can compute 

< T2fc+1 >= 

< r 2i >= jf LgJM iy = jf ( 2 * ) rfgi* (4.2) 

where L = Q + ai(x)C _1 is the genus continuous version of Q(l), and the label ( ) denotes the 
coefficient of £°. In ^3] the resolvent is defined as 



oo 

1 



Therefore in the genus case we have* 



N ^ p k 

k=0 ^ 



<r k >o -A (2k)\ 2k x 



>k ^0 \ 



2k+l 



p k+i kUk + lV. 3 p^ 1 

k=0 * k=0 v ; ^ 

= vT^*) =i( P - Vp^4t) (4.4) 

which is the result of [Oj for the resolvent, provided we set t = xg s and u = — . From this we 
can reconstruct the eigenvalue density 

,<») = ^V^ 5 

In the saddle point method this leads to introducing an auxiliary hyperelliptic Riemann surface 

y 2 = p 2 - it (4.5) 

where all data of the model are encoded. 

As we see from this elementary example, starting from our formulas above we can reconstruct 
all the formulas of the saddle point method. In this case there is no difference between the two 
methods. 



4.2 The cubic model 

The full 7^3,2 model has been discussed at length in [22] and, in particular, in Here we would 
like to consider its decoupling limit c = and single out the cubic potential part, which amounts 
to considering the cubic one-matrix model M.3. In the genus limit this model is described by 
the discrete algebraic equations 

ai = ~2 a °~3^ a ° (4-7) 



*In the continuous limit the couplings get renormalized, t k — ► t& = t k /N and, due to the expansion (|3.15() . we define 



< * >o= tr- 



io 



where, for simplicity and without loss of generality, we have set t\ = 0. One can extract from 
these equations ao and a\ and calculate all the correlators with the algorithm described in the 
previous section. Here we are not interested in this, but rather in analyzing ea. ()4.6|) and its 
solutions. 

In the large TV" limit we shift to x = -8* and, in order to make contact with section 4 of [23] 
for a comparison, we simplify a bit further our notation setting t 2 = — y and £3 = —Ng, where 
g is the cubic coupling constant there. Moreover we denote z = 3gao- Then ea. ()4.6|) becomes 

18g 2 x + z(l + z)(l + 2z) =0 (4.8) 

This can be solved exactly for z and gives the three solutions 

1 1 I(x) , , 

1 l + iy^ l-iV3I(x) , A , 

~~ 2= -2 + ^r + 12 (4 - 10) 

1 l-n/3 l + i\/3J(a;) 11N 

- 3 = -2 + ^r + 12 (411) 

where 

I(x) = 3 1/3 (-324a 2 x + v / 3 V / -l + 34992 5 4 a; 2 ) (4.12) 

From these we can extract three solutions for ao and, consequently, for a\. For small x the three 
solutions can be expanded as follows 

Zl = -18g 2 x - 972aV - 93312a 6 x 3 - 11022480a 8 x 4 + 0(x 5 ) (4.13) 
z 2 = -1 - 18o 2 x + 972o 4 x 2 - 93312o 6 x 3 + 11022480o 8 x 4 + 0(x 5 ) (4.14) 

z 3 = -- + 36o 2 x + 186624o 6 x 3 + 0(x 5 ). (4.15) 

The best way to analyze these solutions is to notice that they represent a plane curve in 
the complex z,x plane. It is a genus Riemann surface with punctures at x = and x = 00, 
made of three sheets joined through cuts running from z = -l/( v / 3108o 2 ) to z = l/(108\/3o 2 ). 
The solutions (|4. 1314. 1414. 15|) correspond to the values z takes near x = 0, away from the cuts. 
In order to pass from one solution to another we have to cross the cuts. We call the Riemann 
surface so constructed the quantum Riemann surface associated to the model}. This Riemann 
surface picture is the clue to understanding the solutions with multiple brane configurations, as 
was explained in [I]. 

Let us analyze the meaning of these three solutions. To this end it is useful to make a 
comparison with [23]: we see that the first solution corresponds to the unique solution found 
there, which corresponds to the minimum of the classical potential (see below). In fact the 
correspondence with |23j can be made more precise: one can easily verify that eqs.(46) there are 
nothing but eas. ()4.6l4.7|) . provided we make the identifications: a + b = oq and (b — a) 2 = Aa\ 
and the rescalings ao — ► v^Ch^i — ► xa\ and g — > g/y/x. In [23] the interval (2a, 26) represents 
a cut in the eigenvalue A plane. Also in [23], as in the previous subsection, we have therefore an 



Un the example of the previous subsection the quantum Riemann surface was just one point. 
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auxiliary Riemann surface. The latter is not to be confused with the quantum Riemann surface 
defined above, although the two are related. 

Let us now discuss the correspondence between our three solutions and the classical extrema 
of the potential. The classical potential for the continuous eigenvalue function \(x) (which is 
\ n /y/N in the large N limit), is V c i = ^A 2 + gA 3 . It has extrema at A = and A = —l/3g. To 
find the classical limit in our quantum approach instead, we rescale tk, k = 2,3 as: i& — ► tk/h, 
and take H — ► 0. This amounts to dropping the last term in ea. (|4.6|) . The extrema are three, 
z = 0,-1,-1/2, which corresponds to ao = 0, — l/3g, — l/6g, not two as in the classical case. 
z = corresponds to the minimum of the potential, z = —1 corresponds to the maximum, while 
z = —1/2 to the flex. The latter solution does not have a classical analog. 

Discussion. Before taking seriously this new quantum solution, we have to check whether 
our method of solving the quantum problem has any flaw. Our method consists in solving the 
quantum equations of motion l|3.8|) . One might wonder whether there are other independent 
conditions beside (|3.8() that one must impose on the solutions. Eas. l|3.8[l are obtained from 
the identity 

J dX 1( iX 2 A (U^)e Vl(Xl)+V2(X2)+cXlX2 vU^ = (4.16) 

and a similar equation for the derivative with respect to A2. Of course one has also 

J d\ l( i\ 2 A (A^r 1 (Ai)e yi ( Al ) +1/2 ( A2 ) +cAlA2 r /m (A 2 )) = (4.17) 

for any k. The question is whether these equations imply additional constraints on the solutions 
of the quantum equations of motion. However one can easily see that eq. I|4.17|l can be written 

as 

Q(l) k (P°(l) + V[{Q{1)) + cQ(2)) = (4.18) 

Therefore a solution to eas. (|3.8|) is automatically a solution to (|4.18|l . Thus these additional 
equations (as well as all those that can be obtained by inserting, instead of A^, generic mono- 
mials of Ai and A2 in (|4. 1 T|t ^1 . cannot further constrain the solutions to (|3.8|) . By the way, 
the W-constraints are obtained precisely by taking the trace of all the expressions like l|4.18|l 
and expressing them in terms of derivative of InZ. Therefore the three quantum solutions we 
have found above are also solutions to the W-constraints (or loop equations). Finally if we 
consider higher order derivatives with respect to Ai in l|4.16[) instead of the first order one, it 
is easy to see that we do not get any additional constraints either. 

This result is somewhat puzzling, but we should remember that the saddle point method is 
semiclassical: one cannot exclude that the quantum problem admits solution without classical 
analog. This is precisely what happens in the present case. One can phrase it also by saying 
that, in general, the large N limit and the h — > limit do not commute. 

The next question is: what is the meaning of the third solution, z = —1/2? Let us recall 
what the other two solutions at z = and z = —1 mean. On the basis of the discussion in 
section 2, we know that they represent two Riemann spheres located at the minimum and the 
maximum of the potential. They replace the continuous family of P 1 which characterizes the 
conifold geometry before the deformation W is introduced. This is the interpretation on the 
basis of classical geometry. What we learn now is that solving the quantum problem we obtain 
a third solution, which we can interpret as a quantum P 1 , located at the flex of the potential. 
This is a pure quantum geometry effect. 
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Before we end this section we would like to make a few remarks. First we notice that the 
classical extrema are characterized by the fact that a\ = 0, while the pure quantum solution 
corresponds to a non- vanishing a\. Moreover, after setting a\ = we get for ao an equation 
that coincides with the classical eigenvalue equation. From this simple example we learn three 
important pieces of informations. 

• The number of solutions of the quantum problem (i.e. the number of solutions to ea. (|4.8(0 
is in general larger than the number of the extrema of the classical potential. 

• The field ao can be regarded as the quantum version of the classical eigenvalue function. 

• The classical extrema are obtained by neglecting the n term in eq. (j4.fi)) and setting at = 0. 

These conclusion are valid in general, except for the fact that, the condition a\ = in the last 
remark must be replaced by the fields a\, a2, b\, 62, being set to zero in the general case. 



4.3 The M 3 , 3 model 

We study the model in the case t\ = si = and limit ourselves to writing down the genus 
quantum equations of motion: 

3i3cao + 2i2cao — 36s 3 t 3 boR + c 2 bo — 12s2t 3 R = 
3S3C&0 + 2s2cbo — 36s3i3do-R — 12s 3 t2R + aoc 2 = 

nc + Rc 2 - l8s 3 t 3 R 2 - 36s 3 t 3 a b R - l2s 2 t 3 a R - 12t 2 s 3 b R - As 2 t 2 R = (4.19) 

01 = —^^-boR — —^R, (12 = — ^^-R 2 

c c c 

6l = _^3 aojR _^ b 2 = -^R 2 

c c c 

From the previous section we see that we are indeed interested in finding all the solutions that 
have an analytic expansion in x = n/N around x = 0. In order to compute all possible solutions 
of this type, that is the quantum vacua, we therefore drop the first term in the lhs of (|4.19|) and 
solve the resulting system. The third equation, in particular, admits the solution R = 0. 

R = (4.20) 
3^0 + 2t 2ao + cb = (4.21) 
3s 3 ^o + 2^260 + ca = (4.22) 

which give rise to four (in general) distinct solutions. The alternative R 7^ leads to 
9 £3^3 a 2 + 3( £3^ c + 8s2t3)ao + 108 S3t3ao6o + 54 s J^l b 2 + i62^a 6^ 

2 S2 S2 S2 S2 

= 3— (c 2 - 16s2t 2 ) + c 2 - As 2 t 2 , (4.23) 

S2 

lc^b 2 + 3(^c + 8s 3 t 2 )a + 108 s 3 Wo + 54^ag + l 6 2^ag6 
z 12 12 12 02 

= 3-(c 2 - 16s 2 t 2 ) + c 2 - 4s 2 t 2 , (4.24) 

This leads to an algebraic equation of order 10 for ao, for instance. Therefore, generically, we 
have 10 (possibly complex) solutions for ao, each of which gives rise to two different values for 
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&o- Altogether we are going to find 24 different quantum vacua. Once again it is interesting 
to compare these solutions with the classical ones. To this end in the above equations we set 
a 2 = b 2 = a i = bi = as well as R = 0, from which we get eas. (|4.2H4.22j) . 

From the first we can get bo = — ^(2>t 3 a^ + 2t 2 a>o), whence we get either ao = or the cubic 
equation 

27s3*i«o + 36t 2 s 3 t 3 al + (12« 3 t| - 6cs 2 t 3 )a + c(c 2 - As 2 t 2 ) = 

Therefore in general we have four classical extrema. In ref.|3j it is shown how to find an explicit 
series expansion in x about each of these solutions. 

Before we end this section it is interesting to discuss the geometric meaning of the first three 
equations in (|4.19|) . We can think of the third equation as a definition of the complex x = n/N 
plane. The two remaining ones are quadratic equations in ao,bo,R. Introducing homogeneous 
coordinates they can be seen to represent two hypersurfaces in P 3 . The intersection is a genus 1 
Riemann surface. 

4.4 The M PljP2 model 

In the general case the matrix rank for Q(l) and Q(2) was given above and the quantum EoMs 
become of course very complicated. It is however simple to write down the equations that identify 
the extrema with classical analog. They are 

V{{oq) + cb = 0, V 2 \b ) + ca = (4.25) 

while all the other fields are set to zero. We have in general (pi — l)(p 2 — 1) solutions of this 
type in perfect correspondence with the classical analysis. The simplest solution is ao = &o = t'- 
other solutions may be hard or even impossible to determine explicitly. Anyhow, once one such 
solution is known it is possible to find explicit expressions for the fields around it in terms of 
x = n/N. The quantum Riemann surface corresponding to this model has the same general 
structure as the one discussed at the end of the previous subsection. It is an intersection curve 
of two hypersurfaces in P 3 . Its genus can be calculated with standard methods in intersection 
theory. For instance, the quantum Riemann surface for the A4^ 3 model is in general genus 2. 



5 Topological open string expansions 

As we have pointed out in section 3 the matrix formalism allows us to compute both closed and 
open string amplitudes. In this final section we would like to show how to obtain the latter. To 
avoid clogging the text with cumbersome formulas we work in the M 2<2 model. The most explicit 
formulas for the latter can be found in [22]. We adapt the results to the present situation by 
rescaling all the coupling, such as tk — > tk/g s , etc. For instance the exact one point function is 

r l 



~\2- k ( N 

(r - 21, 

2 l=0 k=0 

where 



A = ^f T , B= 2 -¥^l (5.2) 
<r - 4s 2 t 2 cr - 4s 2 t 2 
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Now, for large N it makes sense to expand the binomial coefficient in ()5.1|) in powers of N: 




(5.3) 



where 



aw 



E 



1 < k < s 



Pa(a) = 1 



Pk( s ) = otherwise 



l<Sl<S2...<Sfe<S 



As a consequence we get 

r 1 l - k+1 r !9~fc (_i\p+l-k+l 

< t >= E E E (r -ai) IM (,:*)i(i-* + i)i *^'(' - fc) ^ B " 2 ' 9i ^ (5 ' 4) 

2Z=0fc=0 p=0 v ; v ; v ' 

Setting p = h and / = 2g — 2 + h, we can easily extract the contribution to the correlator from 
world-sheets of genus g with h boundaries. This of course requires I — p to be even. But in 
(|5.4[) we have also contributions with odd I — p. We do not have an interpretation for these 
additional contributions. They may be perhaps related to the possibility of describing the effect 
of punctures. 
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